The Structure of the Aoki Phase 
at Weak Coupling 



R. Kenna, C. Pinto and J.C. Sexton, 
School of Mathematics, Trinity College Dublin, Ireland 

January 2001 
Abstract 

A new method to determine the phase diagram of certain lattice fermionic field theories in 
the weakly coupled regime is presented. This method involves a new type of weak coupling 
expansion which is multiplicative rather than additive in nature and allows perturbative calcu- 
lation of partition function zeroes. Application of the method to the single flavour Gross-Neveu 
model gives a phase diagram consistent with the parity symmetry breaking scenario of Aoki and 
provides new quantitative information on the width of the Aoki phase in the weakly coupled 
sector. 



In lattice field theory, there has been considerable discussion on the phase diagrams of theories with 
Wilson fermions (see, e.g., ^ ||, ||, [7|, |8|, ^, 10, These can be considered as statistical 

mechanical systems, and have rich phase structures whose existence is due to lattice artefacts. 
The Wilson fermion hopping parameter is k = 1/2 (Mo + d) where Mo is the dimensionless bare 
fermion mass and d the lattice dimensionality. It is well known that a system of free Wilson 
fermions exhibits a phase transition at k = 1 / (2d) and that massless fermions appear at this point 
in the continuum limit. Discussions concern the extent to which this phase transition persists in 
the presence of a bosonic field. In QCD, where Wilson fermions couple to gauge fields with a 
strength given by the dimensionless coupling g, there are two candidates for the phase diagram. 
In the first, pioneered by Kawamoto, the expectation is that there is a line of phase transitions 
(the "chiral line") extending from the strong coupling limit to the weakly coupled one and along 
which the pion and quark masses vanish |IJ. Such vanishing is symptomatic of spontaneous chiral 
symmetry breaking. Approaching the continuum limit, at g = 0, along the chiral line in particular 
is then expected to recover massless physics. This is still sometimes referred to as the 'conventional' 
picture. 

The second candidate phase diagram for QCD was determined by Aoki on the basis of comparison 
with the Gross-Neveu model 0]. The Gross-Neveu model serves as a prototype for QCD [Oj. 
Indeed, except for confinement, it has features similar to QCD. One of these features is asymptotic 
freedom, so that in the Gross-Neveu model, as in QCD, the continuum limit is taken in the weakly 
coupled zone. Two features distinguish Aoki's phase diagram from the earlier 'conventional' picture. 
Firstly, instead of a single critical line, Aoki's analysis advocates the existence of two lines extending 
from the strongly to weakly coupled limits, with a number of critical points at g = linked by 
cusps (see Fig. 1). The region above the cusps and between the two extended lines is often referred 
to as the Aoki phase. Secondly, the existence of the phase transition in Aoki's scenario is due 
to spontaneous parity symmetry breaking within the Aoki phase, as opposed to chiral symmetry 
breaking. (Indeed, Wilson fermions explicitly violate chiral symmetry.) This is signaled by a non 
zero vacuum expectation value of the pseudoscalar operator ir = tpi'ysifj in the thermodynamic limit. 
The masslessness of the pion is then attributed to the divergence of a correlation length associated 
with this second order phase transition. In the multiflavour case, flavour symmetry is also broken 
in the Aoki phase since the pion, whose expectation value is nonvanishing, also carries flavour. 
The continuum limit has to be approached from outside the Aoki phase since parity and flavour 
are conserved in the strong interaction. The physical meaning of an approach to the continuum 
limit from within the Aoki phase is unclear j2|. There exists substantial evidence supporting Aoki's 
scenario in the strongly coupled regime j| |, f§ |, § 0]. In the weakly coupled regime the evidence 
has, however, been controversial [Q] (see || for recent discussions on this topic). 

In asymptotically free theories the weakly coupled region is the appropriate one for the continuum 
limit. Recently, Creutz |l(J questioned whether the Aoki phase, pinched between the arms of cusps, 
is "squeezed out" at non-zero coupling or whether it only vanishes in the weak coupling limit (see, 
also, 

The purpose of this paper is twofold. We introduce a new type of expansion which is multiplicative 



rather than additive in nature and from which information on the partition function zeroes of the 
theory can be extracted in a rather natural way 0. Secondly, we address the question of the 
"squeezing out" of the Aoki phase at weak coupling. This multiplicative approach to the single 
flavour Gross-Neveu model, shows that the width of the central Aoki cusp is 0(g 2 ) while the Aoki 
phase has not yet emerged at this order from the left and right extremes. 

The Gross-Neveu model is actually a two dimensional model of fermions only, which interact 



through a short range quartic interaction [12|. In Euclidean continuum space, the model with 
a single fermion flavour is given by the four-fermi action 
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where 75- = i *7i72 and the fermion field has 2 spinor components. Bosonizing the action gives for 
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and where <f>(x) and tt(x) are auxiliary boson fields. The corresponding Wilson action in terms of 



dimensionless lattice quantities is S 
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and where the auxiliary fields have been rescaled <fi — ► gcp, ir — > ^7r to explicitly display the order 
of the interactive part of the action. Here, lattice sites are labeled n M = —N/2, . . . , N/2 — 1, where 
N is the number of sites in each of the two directions. We assume N is even. 

Using the lattice Fourier transform, ip(n) = (1/Na) 2 J2k V'(^) ex P {ikna), where a is the lattice 
spacing, the fermionic action can be written 
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where the 27V 2 x 2N 2 fermion matrix is M^(p,q) = M^(p,q) + M^ nt \p,q), with free part 
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and interactive part 
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It is appropriate to impose antiperiodic boundary conditions in the temporal (1-) direction and 
periodic boundary conditions in the spatial (2-) direction in coordinate space. With these mixed 



boundary conditions the momenta for the Fourier transformed fermion fields are p^ = 2-Trp^/Na, 
where p x G {-N/2 + 1/2, . . .,N/2 - 1/2} and p 2 G {-N/2, . . .,N/2 - 1}. Integration over the 
Grassmann variables gives the full partition function 

Z = J T><f>VTrV$V^exp (-^ W) ) oc (detM^) oc (j[X a (p)^ , (9) 

with X a (p) the eigenvalues of the fermion matrix and the expectation values being taken over the 
bosonic fields. 

In the free field case the eigenvalues of are easily calculated and found to be 
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are the Lee- Yang zeroes of the free theory [13]. Note that the eigenvalues (|l0|) and the zeroes ( JTT| ) 



are degenerate with respect to p^ — > — pu- Furthermore, the lowest zeroes in the free case, and those 
responsible for the onset of critical behaviour, are two-fold degenerate in two dimensions. These 
lowest zeroes are 770, (± |pi | , ^2) where |pi| = (N — l)/2 or 1/2 and p2 = —N/2 or and impact on 
the real 1/2k axis at —2, and 2. Finally note that the zeroes in the upper half plane are given 
by a = 1, while their complex conjugates correspond to a = 2. 

The standard additive weak coupling expansion of the full fermion determinant is the Taylor ex- 
pansion of 

detMW = detMW xdet(M(°) _1 M( w )) 
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This expansion is 
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where the indices i and j stand for the combination of Dirac index and momenta (ot,p) which label 
fermionic matrix elements, so that M^ mt ^ represents (\a\p)\M^ mt ^ (p, q)\Xa (q)}. Here \Xp(q)) 
represents a free fermion eigenvalue. The traces in ( |I3| ) are just the diagrams which contribute to 
the vacuum polarization tensor. 

Setting U = (M| nt) ) and % = = (M^ int) M^ nt) ) - (M^ nt) MjJ nt) ), the ratio of partition functions 
is, from (|l~3l), 
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We note at this point that this expansion is analytic in 1/2/c with poles at 1/2k = r{f\ 



For the Gross-Neveu model, calculation of the pure bosonic expectation values is particularly simple. 
Indeed, one has that (4>(n)) = (ir(n)) = and 



(</>(n)0(m)) = (7r(n)7r(m)) = 25 nm . (15) 
The required bosonic expectation values of the matrix elements are found to be 

U = t a<p = , (16) 
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An alternative formulation of the partition function may be obtained by writing the Wilson fermion 
matrix as Af( w ) = 1/2k + H where H is the hopping matrix. The fermion determinant detM^ w ^ = 
det(l/2w + H), is a polynomial in 1/2k since for finite lattice size these matrices are of finite 
dimension. Indeed, for an N x N lattice this polynomial is of degree 2N 2 . Therefore the bosonic 
expectation value of the fermion determinant is also a polynomial of the same degree in 1/2k and 
as such may be written in terms of its 27V 2 zeroes, now labeled r\i. 



We may thus write a 'multiplicative' weak coupling expansion as 
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where Aj = rn — rji^ are the shifts that occur in the zeroes when the bosonic fields are turned on. 
These are the quantities to be determined. Note that the expression © is, like (@), analyt ic m 
1/2k with poles at rjf . Expanding ( Jl8| ) gives 

/detM( w )\ A, l 2 ^ 22 ^ 2 A, A 



M(0 ) > 1 E A (o i ) + 2 EE 

/ 4 = 1 A j 1=1 Jy±l 



dviM^ I ' ^ 2 ^ — A (o) A (°) + -" ' (19) 



In the free fermion theory, the eigenvalues and zeroes of (|T(j) and (11) are two- or four- fold 
degenerate with respect to momentum inversion. Let {n} denote the n th degeneracy class, so that 
the D n eigenvalues Xn{ = . . . = Xn), n are identical to Xn \ say. Let 1/2k = r/i°^ + e and expand 
the additive and multiplicative expressions ( |l4| ) and (|l9|) order by order in e _1 . Indentification of 
the expansions yields relationships between the known quantities ti, and i« and the shifts in the 
positions of the zeroes, Aj, to C(e _1 ) and 0(e~ 2 ). The C(e _1 ) relationship is 

E aJi- E ^io)}= E E -totSt ' ( 2 °) 

having used (16), while that to (D(e~ 2 ) is 
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These relationships can be considered order by order in the coupling as well. Let Aj = r]^ + 77} + 
0(g 3 ), where 77- and and rj^ are, respectively, the order g and order g 2 shifts in the i th zero. One 



finds that the 0(e 1 ) equation to order g is 



while to order g 2 it is 
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Also, the 0(e ) equation, which is entirely 0(g ), is 
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With relations (22)-(|24|), the multiplicative expression ( |l9| ) recovers (14) to 0(g 2 ). Now the additive 
and multiplicative expressions (|l4| ) and (|l8| ) coincide to 0(g 2 ) everywhere in the complex hopping 
parameter plane and arbitrarily close to any pole. 

In the free case, the zeroes responsible for criticality are two fold degenerate. For weak enough 
coupling, one expects these zeroes to govern critical behaviour in the presence of weakly coupled 
bosonic fields too. From (22) and p4]), the first order shifts to two-fold degenerate zeroes are 



T$=±y/i^; , (25) 

where rii € {n} for i = 1 or 2. The second order equation in the two-fold degenerate case is 
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Removing the bosonic field expectation values converts the problem into the determination of the 
eigenvalues of a weakly perturbed matrix whose free eigenvalues are two-fold degenerate. More 
explicitly, with boson expectation values removed, the eigenvalues are Aj = \f* — rj^ — ■q'p , which 
may be determined from conventional time independent perturbation theory. This condition yields 
enough to fully determine the zeroes to order g 2 . Indeed, the second order shifts are 

J 2 ) = V ^ (27) 
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Now using (|T7|), the O(g) and 0(g 2 ) shifts for the erstwhile two-fold degenerate zeroes, rj a (±\pi\,p2) 
(for p2 = or —N/2), are, respectively, 
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The partition function zeroes are 'protocritical points' [14] whose real parts are pseudocritical 
points. In the thermodynamic limit these become the true critical points of the theory and their 
determination amounts to determination of the weakly coupled phase diagram, the critical line 



being traced out by the impact of zeroes on to the real hopping parameter axis. Thus, the phase 
diagram is given to order g 2 by the limit 



i-= lim {v^(p) + V^(p)+V^(p)\ , 
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where p is the momentum corresponding to the lowest zeroes. 



At order g°, the zeroes (11) impact on the real 1/2k axis at —2, and 2, giving three different 
continuum limits, corresponding to the nadirs of the three Aoki cusps (2|. The true continuum 
limit is 1/2k = 2. From (|28| ) and (p9|), the 0(g 2 )-shift is the shift in the average position of the 
two zeroes while their relative separation is O(g). In the thermodynamic limit, the O(g) terms in 



(28) vanish. One finds, numerically, that the imaginary contribution to the 0(g 2 ) term ( |2~9|) also 
vanishes while the real part becomes an iV-independent constant. Indeed, the factor 
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approaches approximately 0.77 and —0.77 for (pi,p2) = (±1/2,0) and (±(N — l)/2, —N/2) respec- 
tively, and corresponding to the rightmost and leftmost critical lines. Also, (|3l| ) is approximately 0.2 
and —0.2 for (pi,P2) = (±{N — l)/2,0) and (±1/2, —N/2) respectively, these two lines generating 
the inner cusp. 

Therefore, the degeneracy of the free fermion critical point corresponding to the central cusp in 
Aoki's phase diagram is indeed lifted and two critical lines emerge in the presence of weak bosonic 
coupling. These critical lines are 1/2k ~ ±0Ag 2 . The Aoki phase does not yet emerge to 0(g 2 ) 
from the left- and rightmost critical points. This is the answer to the question posed by Creutz in 
[jio| at least for the single flavour Gross-Neveu model. Instead the right- and leftmost critical lines 
are 1/2k ~ ±(2 — 1.54<? 2 ) respectively. 




1/(2 K) 

Figure 1: The phase diagram for the Gross-Neveu model in the weakly coupled region (to 0(g 2 )) 
(dark lines) and a schematic representation of the expected Aoki phase diagram (light curves). 



This weakly coupled phase diagram is pictured in Fig. 1 (dark lines). From left to right, the critical 
lines are traced out by the thermodynamic limits of the zeroes indexed by (pi,P2) = (±(A^ — 



l)/2,-iV/2), (±(N - l)/2,0), (±1/2, -N/2) and (±1/2,0) respectively. The lighter curves are a 
schematic representation of the expected full phase diagram. 

In conclusion, we have developed a new type of weak coupling expansion which is multiplicative 
rather than additive in nature and focuses on the Lee- Yang zeroes, or protocritical points, of a 
lattice field theory with Wilson fermions. This expansion is applied to the Gross-Neveu model, 
where the existence of an Aoki phase was first suggested. The weakly coupled regime is the one 
of primary interest as it is there, as with all asymptotically free models, that the continuum limit 
is taken. The method, applied to the single flavour Gross-Neveu model, yields a phase diagram in 
this region which is consistent with that of Aoki and the widths of the Aoki cusps are determined 
to order g 2 . 
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